Bloch-mode extraction from near-field data in periodic waveguides 
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We demonstrate that the spatial profiles of both propagating and evanescent 
Bloch-modes in a periodic structure can be extracted from a single measure- 
ment of electric field at the specified optical wavelength. We develop a system- 
atic extraction procedure by extending the concepts of high-resolution spectral 
methods previously developed for temporal data series to take into account 
the symmetry properties of Bloch-modes simultaneously at all spatial loca- 
tions. We demonstrate the application of our method to a photonic crystal 
waveguide interface and confirm its robustness in the presence of noise. 
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Periodically modulated optical waveguides offer new 
possibilities for controlling the propagation of light. Res- 
onant scattering from periodic modulations can be used 
to tailor the dispersion, enabling in particular a dra- 
matic modification of the group velocity and realization 
of slow-light propagation. Such fundamental effects can 
be directly visualized in experiment with near-field mea- 
surements, which can be used to recover the amplitude, 
phase, and polarization of the electric field at all spa- 
tial locations in the plane of the waveguide [1]. This 
information can then be used to extract the dispersion 
characteristics of the guided modes. 

A commonly used approach to the dispersion extrac- 
tion is through the spatial Fourier-transform (SFT) of 
the field profiles, since peaks in the Fourier spectra cor- 
respond to the wavenumbers of guided modes @, 0, 3- 
However, there exists a fundamental limitation on results 
obtained with SFT: Ak > 2tt/L, where Ak is the resolu- 
tion of the wavenumber, and L is the structure's length. 
Therefore, accurate dispersion results can only be ob- 
tained for long waveguides, extending over many periods 
of the underlying photonic structure. Another limitation 
of the SFT method is that it cannot provide informa- 
tion on the dispersion of evanescent waves, which may 
play an important role close to the structure boundaries 
or interfaces between different waveguides. For example, 
evanescent waves enable efficient excitation of slow-light 
waves without a transition region [Bj]. 

Alternative methods for dispersion extraction have 
been developed to overcome the shortcoming of the SFT 
method. It was shown that an interference of two 
counter-propagating modes can be used to extract their 
wavenumbers [6|], however this technique is not applica- 
ble under the presence of multiple propagating modes or 
evanescent waves. Recently, it was demonstrated that 



dispersion extraction in multi-mode waveguides with in 
principle unbounded resolution is possible even for short 
waveguide sections 0, H|, using approaches based on 
an adaptation of high-resolution spectral methods pre- 
viously developed for the analysis of temporal dynam- 
ics 0, [To| . In this work, we introduce an important gener- 
alization of such methods taking into account the spatial 
symmetry properties of modes in periodic waveguides. 
We show that beyond the dispersion relations, it is pos- 
sible to extract the spatial profiles of all guided modes. 
Our method is applicable to an arbitrary combination 
of propagating and evanescent waves. We illustrate the 
application of this general approach by analysing light 
dynamics at an interface between photonic crystal wave- 
guides designed for coupling into a slow-light mode [5] 



Let us consider a periodic waveguide section, where 
the light propagation in a particular frequency range is 
primarily determined by a finite total number of guided 
modes (M). The value of M can be established based 
on numerical modelling, taking into account both prop- 
agating and evanescent waves. Since each of the modes 
of a periodic waveguide satisfies the Bloch theorem [~Llj , 
the complex electric field envelope of a waveguide mode 
with the index m at the frequency uo can be expressed as 
^m(r; uj) exp(ik m z/d). Here k m are the complex Bloch 
wavenumbers, r = (x, y, z) where x and y are the or- 
thogonal directions transverse to the waveguide and z is 
the direction of periodicity, d is the waveguide period, 
and ifj m are the periodic Bloch- wave envelope functions: 
ipm(z) — ipm(z + d). Then, the total field inside the 
waveguide can be presented as a linear superposition of 
M propagating modes with amplitudes a m and small ad- 
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Fig. 1. (a) Schematic of a two-dimensional photonic crys- 
tal waveguide. The Bloch-wave extraction is performed 
in the PC2 section, (b) Numerically calculated amplitude 
(left) and phase (right) of complex electric field profile at 
the normalized frequency d/\~ 0.2662. 

ditional contributions w(t;uj): 

M 



m=l 

Here w(r;u) can account for the radiation field due to 
the excitation of non-guided waves and for vanishingly 
small evanescent waves which are excluded from consid- 
eration, and this term can also appear due to noise in 
experimental measurements. 

We now present the procedure for the simultaneous 
extraction of the wavenumbers of profiles of the guided 
modes. Let us separate the spatial domain in a number of 
unit cells, (x,y, z + n-d). Here d is the period in the direc- 
tion along the waveguide, zq belongs to a single unit cell 
(z m i n < z < z m i n +d), and n = 1 : N where N is the num- 
ber of cells in the waveguide section. Let us denote with 
U n {r\ u) = E(x, y,z + n-d; u), A m (r; u) = a m ^ m (r; a;), 
and w n (r;uj) = w(r + n • d;u), where r belongs to the 
first unit cell. Then Eq. (pQ) can be written as, 

M 

U n (r; uj) = A m (r; uj) ex.p(ik m n) + w n (r; w), (2) 

m=l 

where we take into account the periodicity of Bloch-wave 
envelopes. If one considers this relation only for a single 
point r in the unit cell, it becomes mathematically equiv- 
alent to the problems considered in spectral analysis of 
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Fig. 2. Complex wavenumbers k m extracted with spatial 
spectral analysis. Propagating and evanescent modes are 
marked by circles and triangles, respectively. Solid curves 
show numerically calculated dispersion. Real and imagi- 
nary parts of each mode are indicated above. Horizontal 
lines at d/X ~ 0.2662 mark the frequency at the inflection 
point. 

temporal series 0, [To| , and high-resolution spectral meth- 
ods can be used to extract the mode wavenumbers 0. 
However, the special property of periodic waveguides is 
that Eq. ([2]) shall be satisfied simultaneously for all spa- 
tial locations r in the unit cell. In order to determine 
the values of k m and A m (r;o;) which describe most ac- 
curately the whole measured field, we employ the least 
squares method and seek a minimum of the functional 



E(t;uj) = ^ a m ijj m (r;uj)ex.p(ik m z/d) + w(r;uj). (1) W = f r 



N 



dl 7 / r En=i l^n| 2 dr, where integra- 



tion is performed over the unit cell. For given wavenum- 
bers, the minimum WA({k m }) = miiiAW is achieved 
when dW/dA m = dW/dA* m = 0. It follows that for each 
point r in a unit cell, the optimal amplitudes satisfy the 
linear matrix equation C H • C • A(r) = C H • J7(r), where 
components of vector A(r) are the optimal amplitude val- 
ues, components of the matrix C are C np = exp(i& p n), 
and vector U(r) components are U n (r) for p = 1 : M 
and n = 1 : N. We can show that W^({fc m }) = 
W A=X = l-J r U H (r)-C- A(r) dr/ J r U H (r) ■ U(r) dr. 
The remaining task is to find the absolute minimum 
W m i n = mmk m Wa (note that, by definition, Wa is real 
and positive), and this can be done numerically, for ex- 
ample by using the 'fminsearch' function in Matlab. 

We apply our method to the two-dimensional [r = 
{x,z)\ photonic crystal waveguide shown in Fig. [U(a), 
where we indicate the unit cell with shading and illus- 
trate the definition of U n . Photonic crystals are viewed 
as effective media for generating slow light because their 
dispersion characteristics can be controlled with design 
parameters [12j. However, coupling light into a slow- 
light waveguide can be a key practical challenge due to 
the field mismatch between the incoming mode with high 
group velocity v g and the slow mode. The photonic crys- 
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Fig. 3. Extracted magnitudes and phases of mode profiles 
A m (x, z) at the slow-light frequency based on the original 
electric field data shown in Fig. ffl(b). 

tal waveguide being analyzed here mediates efficient cou- 
pling into and out of a slow mode without any matching 
region, where an evanescent mode (PC 2) helps match 
the fields of the incoming mode (PC 1) and the slow 
mode (PC 2) without carrying any energy itself [5]. We 
aim to extract the individual modes of PC 2 from a nu- 
merical data of E(x,z), such as shown in Fig. [2(b). In 
PC 2, the dispersion relation contains an inflection point, 
and therefore the total number of modes that primarily 
define the light dynamics is M = 6, accounting for the 
forward and backward slow-light (m = 1,2) and evanes- 
cent modes (m = 3, 4, 5, 6). The fact that the waveguide 
has multiple modes and a relatively short length provides 
good testing environment for our method. 

We note that the spatial spectral analysis can bene- 
fit through the application of additional constraints due 
to the symmetry of modes in lossless dielectric wave- 
guides, wavenumbers of which are related as: &2 = — fci, 
&4 = — &3, fcs = /eg, fee = — &5- Therefore, k\ and ks are 
the independent parameters that fully define the mode 
dispersion, and we take this into account for numeri- 
cal minimization of the functional W . The extracted 
wavenumbers k m are shown in Fig. [2](a) and they are in 
good agreement with numerical calculations. Most im- 
portantly, we simultaneously extract the profiles of the 
corresponding Bloch waves, which amplitudes reflect the 



excitation dynamics. The mode profiles A m (x,z) at the 
slow-light frequency d/X ~ 0.2662 are presented in Fig. [3j 
The profiles clearly show the propagating (m = 1,2) 
and evanescent (m = 3,4,5,6) modes. This provides 
an essential insight into the light dynamics, which can- 
not be inferred directly from the field profiles as shown 
in Fig. [2(b), or from their SFT spectra. 

To test our method in a possible experimental situation 
where noise is present, we add a normal distribution of 
random complex numbers to the electric field U n (x,z) 
with a mean of zero and a standard deviation of 0.05. We 
have successfully recovered the mode dispersion under 
the effect of such perturbations. 

In conclusion, we have presented a general approach for 
the simultaneous extraction of wavenumbers and ampli- 
tudes profiles of Bloch waves based on near-field measure- 
ments in periodic waveguides. We have demonstrated the 
application of this method for the simultaneous charac- 
terization of multiple propagating and evanescent modes, 
and its robustness under the effect of noise. This ap- 
proach may serve to provide essential insight into the 
light dynamics in complex photonic-crystal circuits. 
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